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[57] ABSTRACT 

A filtering method and apparatus in which an input signal is 
subjected to an analog filtering step, then an analog to digital 
conversion step and finally a digital filtering step to produce 
a filtered output signal. The pulse response of the digital 
niter is the mirror image of the pulse response of the analog 
filter and the cascaded filtering steps result in a filtered 
output signal which exhibits little or no group delay or phase 
distortion. The digital filter may be a FIR filter. Also 
disclosed is a phase error correction method whereby a pulse 
is inserted into a circuit or data stream prior to a group delay 
imparting device or function, a pulse response of the group 
delay imparting device or function is measured prior to a 
phase correcting device or function, mirror image filter 
coefficients are calculated using a window function and then 
incorporated into a phase correcting filter arrangement. 

12 Claims, 3 Drawing Sheets 
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GROUP DELAY CORRECTION METHOD 
AND APPARATUS 

FIELD OF THE INVENTION 

This invention relates to signal processing arrangements 
generally and particularly to those signal processing 
arrangements which impart group delay errors to the pro- 
cessed signal. 

BACKGROUND OF THE INVENTION 

A t4 time reversal*' filtering method is a method of filtering 
an input signal twice, once in the forward time direction and 
once in the reverse time direction. The first (forward) filter 
exhibits a first frequency-dependent magnitude response and 
a first frequency-dependent phase response. The frequency 
dependent phase response, or group delay, disadvanta- 
geous^ distorts phase dependent information components 
included within the input signal (e.g.. audio or video infor- 
mation included within a television signal). To compensate 
for the distorting effect of the first (forward) filtering, a 
second (reverse) filter is cascaded with the first (forward) 
filter. The second (reverse) filter exhibits a second 
frequency-dependent magnitude response and a second 
frequency-dependent phase response. If the group delays of 
the two filters are complimentary, then the phase-distorting 
effect of the first filter is canceled by a "phase -restoring" 
effect of the second filter. This time reversal filter response 
advantageously provides rninimurn phase distortion of the 
filtered signal. In addition, the magnitude response of the 
time reversal filtered signal is the product of the magnitude 
of one of the forward and reverse filters. 

A time -reversal filter is relatively easy to implement using 
a computer running the appropriate software. For example, 
the MATLAB™ computer program, including the Signal 
Processing Toolbox, from The Mathworks, Inc., includes a 
command **filtfilt w which will perform the above described 
forward and reverse filtering of a signal stored in a com- 
puter's memory. Unfortunately, if it is necessary to perform 
time-reversal filtering of a high-bandwidth signal in real- 
time then a dedicated hardware implementation becomes 
necessary. 

An example of a prior art hardware implementation of a 
time-reversal filter is shown in FIG. 4 and includes an analog 
low pass filter (LPF). an analog to digital (A/D) converter, 
an infinite impulse response (IIR) filter, four last in first out 
(LIFO) memory arrangements and four switches 
(SW1-SW4). An input signal IN is filtered and digitized by 
the cascaded LPF and A/D converter, resulting in a 'for- 
ward" filtered input signal. The forward filtered signal must 
now be Reverse" filtered by the IIR. The **time reversal" is 
accomplished by storing the forward filtered signal in the 
UFO memory arrangement. When the forward filtered sig- 
nal is read form the LIFO memory arrangement the more 
recent samples will be retrieved first, thus producing a 'time 
reversed*' forward filtered signal 

Unfortunately, since the LIFOs cannot hold an infinite 
amount of data the signal stream needs to be partitioned and 
reverse filtered by the IIR segment by segment. The seg- 
menting is accomplished by using the four LIFOs as follows. 
In a first mode of operation all the switches are in the "1" 
position (as shown in FIG. 4). the first and third LIFOs are 
receiving and storing data while the second and fourth 
LIFOs are transmitting data previously stored (i.e.. during a 
second mode of operation). In the second mode of operation 
all the switches are in the "0" position, the first and third 
LIFOs are transmitting data previously stored (i.e.. during 
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the first mode of operation) while the second and fourth 
LIFOs are receiving and storing data. A controller (not 
shown) controls the switches and LIFOs. thus causing the 
modes of operation to be alternated. 

5 Disadvantageous^, as a result of the segmented 
processing, the input signal stream is interrupted and signal 
corruption occurs at the segment boundaries. In a digital 
television system the segment length may be selected to be 
one horizontal line, thereby causing any boundary corrup- 

io tion to occur during the horizontal blanking interval. 
However, if the data stream to be filtered is continuous, such 
as an MPEG data stream (which has no blanking interval) or 
a communication channel, then alternate methods must be 
used. For example, an "overlap-add" technique may be used 

15 to help reduce the corruption. This technique requires con- 
siderable amount of hardware to store data beyond a seg- 
ment boundary and re-integrate the data stream at the exact 
boundary segment. Moreover, the results of this technique 
are not always acceptable. 

20 It is therefore desirable to implement a real-time filtering 
system which provides significantly reduced phase 
distortion, does not require segmentation of an input signal 
and does not require an excessive amount of hardware to 

^ implement. 

SUMMARY OF THE INVENTION 

The subject invention concerns a filtering method and 
apparatus in which an input signal is subjected to an analog 

30 filtering step, followed by an analog to digital conversion 
step and finally a digital filtering step to produce a filtered 
output signal. The pulse response of the digital filter is the 
mirror image of the pulse response of the analog filter and 
the cascaded filtering steps result in a filtered output signal 

35 which exhibits little or no group delay or phase distortion. 
The digital filter may be a FIR filter. The characteristics of 
the digital filter may be periodically updated by exciting the 
analog filter, measuring its pulse response and calculating 
new filter coefficients which, when utilized by the digital 

40 filter, cause the pulse response of the digital filter to be a 
mirror image of the pulse response of the analog filter. 

BRIEF DESCRIPTION OF THE DRAWING 

The invention will be described with reference to the 
45 accompanying Drawing, in which: 

FIG. 1 a block diagram of a quadrature phase- shift keying 
(QPSK) demodulator including a filtering arrangement 
according to the invention; 
^ FIG. 2 is a block diagram of a filtering arrangement 
according to the invention and suitable for use in the QPSK 
demodulator of FIG. 1; 

FIG. 3 is a graphical representation of several waveforms 
useful for understanding the invention; and 
55 FIG. 4 is a block diagram of a prior art time reversal filter. 

DETAILED DESCRIPTION 

FIG. 1 shows a block diagram of a "direct conversion'* 
tuner in which a quadrature phase-shift keying (QPSK) 

60 demodulator is used to convert a radio frequency signal into 
in-phase (I) and quadrature-phase (Q) signals. The direct 
conversion tuner will be more fully described below. This 
type of receiver requires extremely sharp low pass filters 
(LPFs) after the quadrature mixers in order to suppress 

65 adjacent channels. Such filters tend to introduce consider- 
able group delay distortion. Equalization of the phase- 
distorting effects of the LPFs is often critical to maintaining 
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data integrity of the modulated signal. Such phase 
equalization, especially of a signal which cannot be readily 
segmented, is difficult to achieve using time reversal filters. 
As such, a new filter methodology suitable for filtering 
non-segmented signals without introducing phase distortion, 
phase is desirable. Such a filter is shown in FIG. 2 and 
incorporated into the QPSK receiver of FIG. 1. 

The invention is directed toward the problem of group 
delay errors imparted to an information bearing signal or 
data stream by a group delay imparting device or function 
(Le.. error source) having a characteristic phase response. 
The characteristic phase response of the group delay impart- 
ing device or function may be determined by inserting a 
pulse into a circuit or data stream prior to the group delay 
imparting device or function. The pulse response of the 
group delay imparting device or function is measured prior 
to a phase correcting device or functioo. mirror image filter 
coefficients are calculated using a window function and then 
incorporated into a phase correcting filter arrangement The 
group delay errors in the information bearing signal or data 
stream may be substantially reduced by subjecting the signal 
or data stream to a processing step having the mirror image 
phase response. 

Referring to FIG. 2, a block diagram of a filtering arrange- 
ment according to the invention is shown. An input signal IN 
includes a carrier component and an information component 
modulated upon the carrier component. The input signal is 
received by an analog filter, illustratively a low pass filter 
having a characteristic gain response and a characteristic 
phase response. The analog filtered signal at the output of the 
LPF is coupled to a conventional analog to digital (A/D) 
converter. The A/D converter converts the analog-filtered 
signal into a data stream which is coupled to a phase 
corrector. Ideally, the phase response of phase corrector is 
the exact opposite (i.e., mirror image) of the phase response 
of the LPF. Such an ideal mirror image response has the 
desirable effect of exactly canceling any group delay 
imparted to the input signal by the analog filter. 
Advantageously, this mirror-image filtering method does not 
require segmentation of the input signal thereby eliminating 
the boundary problem exhibited by the time reversal 
arrangement of FIG. 4. Note that waveforms C and D of 
FIG. 3 are graphical representations of, respectively, a 
normalized pulse response of an analog low pass filter 
(truncated) and a normalized mirror-image of the pulse 
response. 

While a substantially ideal phase corrector may be con- 
structed using an infinite impulse response (HR) filter, such 
a filter tends to be complex and expensive. A finite impulse 
response (FIR) may also be used; however, some compro- 
mises must be made since an infinite pulse response from an 
FIR requires an infinite (or very large) number of taps. The 
inventor recognized that a reasonable compromise between 
the number of taps and the performance of the FIR phase 
corrector may be achieved by using a method which will 
now be described. 

For the purpose of this discussion it is assumed that the 
analog filter of the filtering arrangement of FIG. 2 is an 
extremely sharp LPF. As an example, assume the LPF is a 
seven pole elliptical LC filter having a -3 dB cutoff fre- 
quency of 2.86 MHz. passband ripple of 3dB. a transition 
bandwidth of 440 KHz and a stop-band attenuation of 57 dB 
beginning at approximately 3 J MHz. The elliptical filter 
also exhibits a group delay nonlinearity having a peak 
nonlinearity at the cutoff frequency of 2.86 MHz, typical of 
elliptical filters. This very sharp LPF is relatively inexpen- 
sive and has a characteristic gain response which makes the 
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LPF suitable for use in the direct conversion receiver of FIG. 
1. Unfortunately, the characteristic phase response of this 
filter causes group delay errors in the information compo- 
nent of the received signal IN. 

5 FIG. 3 is a graphical representation of five waveforms (A, 
B. G D and E) useful for understanding the invention. The 
horizontal axis corresponds to time and the vertical axis 
corresponds to amplitude. To clearly present the five wave- 
forms in a single figure the third (C), fourth (D) and fifth (E) 

io waveforms have been vertically offset from the first two 
waveforms by, respectively, -0.5 units, -1.0 units and -1.5 
units. 

The second waveform (B) of FIG. 3 is a graphical 
representation of a normalized time-domain pulse response 

15 of the analog low pass filter described above with respect to 
FIG. 2. As previously described, the horizontal axis is a 
measure of time and the vertical axis is a measure of 
amplitude variations of a filtered output signal resulting 
from an input pulse signal The pulse response of the ideal 

20 phase corrector for the LPF of FIG. 2 is a mirror image of 
the pulse response of second waveform (B) at the vertical 
axis. Such a phase corrector may be realized by using an FIR ; 
having 140 taps. However, since a 140 tap FIR filter is quite ; 
complex it is desirable to reduce the number of taps by. e.g.. 

25 windowing and truncation. 

Hie first waveform (A) of FIG. 3 is a graphical represen- 
tation of a normalized smoothing window, illustratively a \ 
modified Kaiser window, suitable for use in a filtering 

^ arrangement according to the invention. The modifications 
may be described graphically as truncating the left half of a 
normal Kaiser window and horizontally displacing the cen- 
ter of the remaining right half. The window interval includes 
most but not all, of the pulse response information from the 

35 LPF (see waveform B). 

The third waveform (C) of FIG. 3 is a graphical repre- \ 
sentation of the product of the first waveform (A) multiplied 
by the second waveform (B). Thus, the pulse response 
information from the analog fiber which is within the ; 

40 window (i.e., time 0 up to time 33) is used to calculate the 
filter coefficients for the FIR phase corrector filter. The pulse 
response information from the analog filter which is outside 
of the window is multiplied by zero, Le.. truncated. In the 
exemplary embodiment this calculation is performed in the 

45 digital domain by. e.g.. a mathematics program (e.g., 
MATLAB™) operating in a microprocessor, microcontrol- 
ler or other computing device (e.g., digital signal processor). . 
The third waveform (C) clearly shows a substantial reduc- 
tion in the residual amplitude components (i.e.. ringing) due ; 

50 to the windowing function. This reduction is superior to that 
which may be achieved using simple truncation (i.e.. a 
rectangular window function) since simple truncation would 
leave higher amplitude components near the zero 
multiplication, or truncation point 

55 The fourth waveform (D) of FIG. 3 is a graphical repre- : 
sentation of the inirror image of the third waveform (C). 
Note that the mirror image of the "truncated" portion of the 
third waveform (D) is not shown. The above-described ; 
windowing method allows for greater truncation of the 

60 number of taps used by the FIR filter of the phase corrector. 
While acceptable performance utilizing only 32 taps is made 
possible by this method, the invention may be practiced 
using a greater or lesser number of taps. The discrete 
component of the fourth waveform (D) of FIG. 4 depicts the ; 

65 normalized pulse response of a 32-tap FIR filter suitable for 
use as a phase corrector for the above-described analog LPF. 
It was experimentally determined that a 40 up filter would : 
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provide very little improvement over the 32 tap filter. this filter adjustment process the performance of the filtering 

However, if less than approximately 28 taps are used then arrangement may be maintained over time and temperature, 

the performance of the phase corrector rapidly drops. For example, the process may be performed each time the 

The fifth waveform (E) of FIG. 3 is a graphical represen- « powered up. after a specific change in temperature has 

tation of a normalized pulse response of the above-described * been measured by the controller <e.g.. ±15° Celsius) or after 

analog low pass filter cascaded with the above-described a certain amount of time has elapsed. 

FIR phase corrector filter. The overall pulse response is Referring now to FIG. 1, the above-described filtering 

fairly symmetrical, though the effect of the truncation can method will now be discussed in the context of a QPSK 

clearly be seen. Note that the right side of the response plot receiver operating according to the direct conversion prin- 

approaches zero more slowly than the left side of the 10 ciple. That is, the QPSK receiver of FIG. 1 demodulates a 

response plot. A filter having an almost symmetrical pulse signal from the RF domain into baseband I and Q signals, 

response has a low group delay or phase distortion charac- There are no image responses because the frequency of the 

teristic. In this case the overall response results in an overall conversion stage is located within the spectrum of the 

phase response which is substantially flat up to about 2.7 desired RF signal. In addition, the very low frequency range 

MHz. which includes 95% of the analog filter passband. 1 5 of the signal produced at the output of the conversion stages 

Thus, the filtering arrangement produces an output signal makes it possible to readily provide sharp filters which can 

OUT in which the phase distortion or group delay errors readily reject adjacent signal in the I and Q channels, 

imparted to the input signal IN by the LPF have been Unfortunately, such sharp LPFs introduce a large amount of 

substantially corrected by the phase corrector. group delay distortions which increase, e.g.. the bit error rate 

The filter coefficients may be calculated using a method 20 « symbol error rate of a received signal QPSK signal. The 

such as the following. First, a portion of the data stream, e.g., &<**P distortions may make proper decoding of the 

the discrete data points of waveform B, is received by the received signal by subsequent processing stages impossible, 

controller and stored in a memory arrangement as a sample In the QPSK receiver of FIG. 1, low pass filters LPF A and 

data set having a plurality of data elements corresponding to LPF B , A/D converters ADC A and ADC B, phase correctors 

discrete data samples. The stored sample data set is then 25 PHASE CORR. A and PHASE CORR. B, resistors RA and 

multiplied by a smoothing function, such as the modified RB and a controller CONTROLLER operate in substantially 

Kaiser function depicted in waveform A. The product of this the same manner as described above with respect to FIG. 2. 

multiplication is a data set having a plurality of data ele- Thus, the overall group delay distortions are canceled in 

ments which correspond to the discrete data samples of substantially the same manner as described above with 

waveform C. The product data set is then truncated to, e.g., 30 respect to FIG. 2. 

32 data elements. The positions of the data elements within Two quadrature local oscillator signals expressed as sintOo 
the data set are then mirror imaged That is. the first and and cosco^ and are applied to, respectively, mixers MA and 
32nd data elements are switched, the second and 31st data MB. Mixer MA produces a channel A intermediate fre- 
elements are switched, etc. The resulting mirror image data quency (IF) signal IF A while mixer MB produces a channel 
set is now used to derive the filter coefficients for the FIR 35 B IF signal IF B. Each of the IF signals IF A and IF B 
phase correction filter. includes a respective double frequency term (Le„ 2o\ > ) as a 
The filter coefficients of the FIR phase corrector filter are consequence of the demodulation process. The purpose of 
determined by the type of analog filter used (i.e.. high pass, the low pass filters LPF A and LPF B is to eliminate these 
low pass, band pass etc.) and the inherent stability of the ^ double frequency terms and, as previously discussed with 
analog filter. For example, in the case where an exception- respect to FIG. 2. to suppress adjacent channels. The pur- 
ally stable filter is used it may be necessary to calculate or pose of the phase correctors PHASE CORR. A and PHASE 
empirically determine the phase response of the filter only CORR. B is to reduce the group delay distortion introduced 
once. However, the components used in the analog filter are into the remaining (Le.. pass-band) frequency terms, 
usually subject to drift and aging problems which are not 45 It will be apparent to those skilled in the art, that although 
easily compensated for. For example, capacitors may change the invention has been described in terms of specific 
value as they age or as the ambient operating temperature examples, modifications and changes may be made to the 
changes. As such, the characteristics of the filters, including disclosed embodiments without departing from the essence 
any phase distortion, will change with, e.g.. time and tern- of the invention. 

perature. ^ For example, the invention has been described in terms of 

To properly compensate for changes in filter characteris- a phase corrector which operatively compensates for group 

tics the phase corrector must also be adjusted. This adjust- delay error imparted to a signal by an analog low pass filter, 

ment is made by periodically measuring the pulse response The error could also have been caused by a digital low pass 

of the analog filter and recalculating the filter coefficients filter, an analog high pass filter, etc It must be noted that the 

used by the phase corrector. The filter adjustment process 55 shape of the window may be adapted for use with other types 

will now be described using the previous LPF filter as an of filters (e.g., more or less sharp LPFs. high pass filters, 

example. band pass filters, etc.). The window function selected was a 

Referring to FIG. 2. a unity pulse is generated by a Kaiser window although there are a number of window 

controller and coupled to the analog LPF by a resistor R. The functions which will provide suitable results. Moreover, 

unity pulse may be coupled to the LPF when. e.g.. the input 60 depending upon the relative complexity of the group delay 

signal IN is either not present or rendered non-interfering by error imparted to the signal by the error source, the number 

switching or other means (not shown). The analog filtered of *ap s and, therefore, filter coefficients required by the 

signal at the output of the LPF is digitized by the A/D phase corrector may be more or less than 32. 

converter (see FIG. 3, waveform B). The digital data stream Finally, it must be understood that the invention is 

is coupled to the controller, which then calculates the filter 65 directed to the broader problem of group delay errors, 

coefficients. The calculated filter coefficients are then regardless of the source of the error. The invention may be 

coupled to the phase corrector. By periodically performing practiced by inserting a pulse into a circuit prior to a group 
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delay imparting device, measuring a pulse response prior to 
a phase correcting device, calculating phase correcting filter 
coefficients using a window function and incorporating the 
coefficients into a phase correcting filter arrangement. 
Moreover, a signal other than a pulse may be used to 5 
produce a response characteristic of the group delay impart- 
ing device. 

These and other modifications are intended to be within 
the scope of the present invention defined in the following 
claims. 10 

What is claimed is: 

1. A signal processing method comprising the steps of: 
receiving an analog input signal; 

processing said analog input signal using an analog filter J5 
to produce a filtered signal, said analog filter having an 
associated phase characteristic; and 

converting said filtered signal to a digital data stream 
using an analog to digital converter; and 

processing said digital data stream using a digital filter to 20 
produce an output signal, said digital filter having an 
associated phase characteristic which is substantially 
the rnirror image of the phase characteristic associated 
with said analog filter. 

2. The signal processing method of claim 1 further 25 
comprising the steps of: 

injecting a pulse signal into an input of said first process- 
ing means; 

measuring the response of said first processing means to 

said pulse signal; 
calculating a mirror image of said response of said first 

processing means to said pulse signal; 
calculating finite impulse response filter coefficients using 

a portion of said calculated rnirror image; 35 
utilizing said calculated coefficients in said digital filter, 

said digital filter being a finite impulse response filter. 

3. The signal processing method of claim 2 wherein: 
said filter coefficients are calculated using a smoothing 

function; and 40 
said portion of said calculated mirror Image is within the 
right half of said smoothing function. 

4. The signal processing method of claim 3 wherein: 

said smoothing function is a modified Kaiser window. 45 

5. A phase error correcting method comprising the steps 
of: 

receiving an input signal having a group delay error 
caused by a first signal processing means, said first 
processing means having an associated phase charao 50 
t eristic; 

processing said input signal using a digital filter to pro- 
duce an output signal, said digital filter having an 
associated phase characteristic which is substantially 
the mirror image of the phase characteristic associated 55 
with said first processing means; 

injecting a pulse signal into an input of said first signal 
processing means; 

measuring the response of said first signal processing 
means to said pulse signal; 



calculating a mirror image of said response of said first 
signal processing means to said pulse signal; and 

adjusting said digital filter in response to said calculated 
mirror image. 

6. Hie phase error correcting method of claim 5 wherein 
said adjusting of said digital filter comprises the steps of: 

calculating finite impulse response filter coefficients using 
a portion of said calculated mirror image; 

utilizing said calculated coefficients in said digital filter, 
said digital filter being a finite impulse response filter. 

7. The phase error correcting method of claim 6 wherein: 
said filter coefficients are calculated using a smoothing 

function; and 

said portion of said calculated mirror image is within the 
right half of said smoothing function. 

8. The phase error correcting method of claim 7 wherein: 
said smoothing function is a modified Kaiser window. 

9. Phase correction apparatus, comprising: 

means for receiving an information signal having a group 
delay error caused by a first signal processing means, 
said first signal processing means having an associated 
phase characteristic; 

digital filter means for filtering said information signal to 
produce an output signal, said digital filter means 
having a phase characteristic which is a substantially 
mirror image of said phase characteristic of said first 
signal processing means, said output signal produced 
by said digital filter means exhibiting a substantially 
reduced group delay error; and wherein 

said first signal processing means comprises analog fil- 
tering means for filtering an input signal; 

said digital filtering means comprises analog to digital 
converter means coupled to said analog filtering means 
for converting an analog filtered input signal to a digital 
data stream. 

10. The phase correcting apparatus of claim 9. wherein: 
said digital filter is a finite input response filter. 

11. The phase correcting apparatus of claim 10. further 
comprising: 

means fox injecting a pulse signal into an input of said first 
signal processing means; 

means for measuring the response of said first signal 
processing means to said pulse signal; 

means coupled to said measuring means for calculating a 
mirror image of said response of said first signal 
processing means to said pulse signal and for calculat- 
ing finite impulse response filter coefficients using a 
portion of said calculated mirror image; and 

means for providing said calculated coefficients to said 
digital filter. 

12. The phase correcting apparatus of claim 11, wherein: 
said calculating means utilizes a smoothing function to 

produce said filter coefficients. 
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